clc;
clear;

plot_basis_function_2d("Q2", 1);

hold on
[X, Y] = meshgrid(linspace(-1, 1, 100), linspace(-1, 1, 100));
z = @(x,y) zeros(size(x,1), size(x,2));
f = z(X,Y);
surf(X, Y, f, "FaceColor", [1 0 0]);
colorbar;